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ABSTRACT 

M-dwarfs in extremely wide binary systems are very rare, and may thus have different formation processes 
from those found as single stars or close binaries in the field. In this paper we search for close companions 
to a new sample of 36 extremely wide M-dwarf binaries, covering a spectral type range of Ml to M5 and a 
separation range of 600 - 6500 AU. We discover 10 new triple systems and one new quadruple system. We 
carefully account for selection effects including proper motion, magnitude limits, the detection of close binaries 
in the SDSS, and other sample biases. The bias-corrected total high-order-multiple fraction is 45^j^% and the 
bias-corrected incidence of quadruple systems is <5%, both statistically compatible with that found for the 
more common close M-dwarf multiple systems. Almost all the detected companions have similar masses to 
their primaries, although two very low mass companions, including a candidate brown dwarf, are found at 
relatively large separations. We find that the close-binary separation distribution is strongly peaked towards 
<30AU separations. There is marginally significant evidence for a change in high-order M-dwarf multiplicity 
with binding energy and total mass. We also find 2cr evidence of an unexpected increased high-order-multiple 
fraction for the widest targets in our survey, with a high-order-multiple fraction of 21^'''% for systems with 
separations up to 2000AU, compared to 11^^^^% for systems with separations >4000AU. These results suggest 
that the very widest M-dwarf binary systems need higher masses to form or to survive. 



1. INTRODUCTION 

Multiple star systems offer a powerful way to constrain the 
processes of star formation. The distributions of companion 
masses, orbital radii and thus binding energies provide impor- 
tant clues to the systems' formation processes. A commonly 
invoked model for very low mass star formation is fragmenta- 
tion of the initial molecular cloud core, followed by ejection 
of the low mass stellar embryos b efore mass ac cretion has 
completed - the ejection hvpothesis dR ei p urth & Clarke 200 1 , 
IClose et al.ll200llSiegler et al.ll2005[ |Konopackv et al. 2007). 
The model successfully predicts the observed field close bi- 
nary frequency and median separation (e.g. Bate 2009), and 
predicts that essentially no binaries with low binding energies 
are expected to survive the ejection. 

Nevertheless, several extremely wide M-dwarf b i naries 
have been detected in the field dPh an-Bao et al.l 120061 
Artigauetal. 2007. ICaballerol 120071" Artigau et al. 2007^ 
Law et all 2008, Radig an et all 120091 iZhang et all |2010[ 
Faherty et a l. 2010). With total system masses of around 



0.2 M0 and orbital radii of thousands of AU, these systems 
have a very low binding energy compared to the limits set by 
the ejection hypothesis. 

On the other hand, there is a sharp binding energ y cutoff 
observed for almost all field VLM binaries (e.g. Clos e et aTl 
l2003L lKonopack y et al.l2007l) . so perhaps the ejection hypoth- 
esis is not invalidated by the few known very low binding en- 
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ergy systems. It is plausible that these "uncommon" systems 
represent a differ ent mode of wide binary production, such as 
mutual ejection (iKouwenhoven et ani201 0). If these unusu- 
ally wide systems are indeed formed in a different way, the 
higher-order multiplicity statistics of the systems may be dif- 
ferent from that of systems with larger binding energies. 

In addition to binding energy considerations for VLM stars, 
empirical limits on binary separation have also been found 
among higher-mass M-dwarfs: log amax — 3 . 3Mfn i +1.1 for 
total masses greater than 0.3Mo (Reidetal. 2001), and for 
lower masses amax - 1400M,^qj (Burgasseret al. 2003), where 
Mtot is the total system mass in solar masses and a^ax is the 
maximum separation in AU. Systems that violate these limits 
could be considered "unusual", and may similarly have differ- 
ent multiplicity properties from the much more common close 
M-dwarf binaries. 

The total masses of such unusually wide systems are gen- 
erally inferred from seeing-limited photometric or spectro- 
scopic fits to determine the spectral types of the bodies mak- 
ing up the binary. The discovery of closer companions can 
thus dramatically alter the inferred total masses of the sys- 
tems, and therefore the expected applicable separation lim- 
its. If a significant fraction of apparently wide binaries are 
actually much higher mass high-order multiple systems, the 
inferences on their formation environments must be corre- 
spondingly altered. In addition, wide binary systems which 
survive the formation process may have experienced smaller 
dynamic influences on their evolution. It is conceivable that 
lower-mass substellar or even planetary companions may be 
more common in such systems. 

New proper m otion surveys bas ed on the Sloan Digital Sky 
Survey (SDSS; I York et al.l l2000h are finding large samples 
of M-dwarf bin aries which app arently violate the Reid et al] 
(2001) & Burgasser et al. ("2005) empirical limits (e.g. SLoW- 
PoKES: Dhital et al. 2010) . For the first time, these new sam- 
ples allow us to undertake detailed statistical studies of the 
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Fig. 1 . — The separations and masses of our wide bina ry sample. Masses are 
inte rpolated from the spectral type / mass relations in 'Kraus & Hillenbrand} 
1200 7a. The empirical field binaiy separation cutoif from Reid et al. (2001) is 
shown as a solid line; we define a system as "unusually wide" if it falls below 
this line. 

systems' multiplicity properties. 

In this paper we target 36 extremely wide binaries, cho- 
sen to cov er the mass and separation range outlined by the 
iReid et all (.200 1) limits. Using Kec k Laser Guide Star 
Adaptive Optics (LGS-AO: IWizinowich et al. 2006) and Palo- 
mar natural guide star adaptive optics (I Trov et al.ii200 0). we 
searched each target system for companions at separations 
ranging from 4 to 100 AU, and at contrast ratios sufficient 
to detect brown dwarf companions. 

We detail the construction of the new sample of extremely 
wide M-dwarf binaries in Section |2] The observations and 
data reductions are described in Section [3] Sections |4] and 
|5] detail and discuss the survey results, and we conclude in 
Section|6l 

2. CONSTRUCTION OF THE WIDE M-DWARF BINARY SAMPLE 

The targets were selected from a preliminary version of the 
SLoWPoKES catalog, a sample of wide (> 500 AU), low- 
mass (mid-K-mid-M) common proper motio n (CPM) binary 
systems in the SDSS Data Release 7 (DR7: lAbazajian et al.l 
12009) . We briefly describe the se lection algorithm h ere; the 
full selection process is detailed in lDhital et al.l (120101) . 

The SDSS DR7 photometric catalog has more than 180 mil- 
lion stellar sources, of which -109 million are low-mass mid- 
K-late-M dwarfs; many sources also have measured proper 
motions in the SDSS/USNO-B matched catalog (Munn et al. 
2004, 2008). To identify CPM pairs from this large database 
Dhital et al.. (.2010) used a relatively high proper motion sam- 



ple (// > 40 mas yr"') to avoid high field contamination. 
Candidate binaries were identified at angular separations of 
7-180", very conservatively requiring photometric distances 
and component proper motions to individually match within 
1-cr [0.3-0.4 magnitudes (.Bochanski et al.. .2010 ) in photo- 
met ric distance modulus, and ~2.5-5 mas yr"' in proper mo- 
tion (iMunn et al.l2004 ','2008')l. A detailed galactic model was 
used to quantify the probability of a random alignment, and 
only pairs with a < 5% chance alignment probability were 
accepted. The final SLoWPoKES catalog contains 1336 very 
wide, low-mass CPM pairs. 

To generate our close-companion-search target list from the 
SLoWPoKES catalog we required 1) both components to be 
M-dwarfs and later than Ml ; 2) the components to have an an- 
gular separation <40"(to allow simultaneous imaging of both 



components); and 3) at least one component to be brighter 
than rx!l7. The latter requirement was imposed due to poor 
conditions during all our LGS AO observations, and effec- 
tively limited the primaries of our targets to be earlier than 
M5. The resulting 36-binary sample (Table [TJ covers a total 
range of 0.35 - 0.85 Mq in mass and 600 - 6500 AU in sepa- 
ration (Figure |2]k 

Using the Dhital et al.l (120101) galactic model and the galac- 
tic positions of our targets, we estimate that our target binaries 
each have a 96% chance of being members of the thin disk, a 
^4% chance of being part of the thick disk, and a < 1 % chance 
of being halo members. A reduced proper motion diagram 
supports this conclusion (Figure|3ll, although the spread in es- 
timated tangential velocities suggests that a few of the targets 
may be members of the thick disk^ 

Our sample is bisected by the IReid et al.1 (1200 Ih field M- 
dwarf maximum separation distribution (Figure[T|l . 23 of the 
36 tar gets are considered unusually wide by the IReid et akl 
(12001b definition. We note that th ere is no clear re duc- 
tion in our target numbers beyond the iReid et al] (1200 Ih cut- 
off", suggesting that the new binary surveys n ow available 
may require the revisiting of this criterion (e.g. ' Dhital et"an 
(2010)). None of our targets violate the suggested alter- 
nate Jeans-length based cri t eria on the system binding energy 
dZuckerman & Son d 120091; iFahertv et aLHIoTO ), although all 
are close to the minimum allowed binding energies and allow 
us to search for a variation in multiplicity properties across a 
wide range of masses and separations. 

3. OBSERVATIONS AND DATA REDUCTION 

We observed the targets with the Keck II and Palomar AO 
systems over 4 nights in 2008 and 2009. Keck LGS-AO ob- 
servations, on 21 May 2008, 8 August 2008, and 13 April 
2009, were undertaken in poor conditions with variable cloud 
cover on all nights. This limited us to the brightest targets 
of our sample, but did not significantly affect the contrast- 
performance-limited minimum detectable companion masses. 
Palomar NGS-AO observations were performed on 8 July 
2009 in exceptional seeing, allowing us to use targets as faint 
as r= 16.5 as natural guide stars. 

We obtained /Tj-band observations at a depth suitable to de- 
tect low-mass brown dwarfs around each target. Stars with 
binary components detected by eye during observations were 
also observed in J and H for photometric characterization if 
conditions and adaptive optics performance allowed. 

The field mid-M-dwarf median binary separation is 
<10 AU (Law et al. 2008; Bergfors et al. 2010). At our tar- 
gets' typical distances of 100 pc, resolving this separation re- 
quires the 40 milliarcsecond angular resolution provided by 
Keek's NIRC2 camera operating in //-band with AO. During 
the majority of the Keck observations conditions were poor, 
only allowing moderate-Strehl Ks observations, and limiting 
the minimum detectable separation to ~8 AU. A few hours 
of the 8 August 2008 night were, however, good enough to 
push to diffraction-limited resolution in //; this subset is large 
enough to allow us to search for any major systematic effects 
stemming from a lack of sensitivity to companions between 4 
and 8 AU. 

The target lists, object properties, image resolutions, and 
observation epochs are summarized in Table [T] Conditions 
were too poor during the observations of SLW1530-(-3935B 
to provide useful limits on its close companion population. 
However, the primary was well imaged and resolved to be a 
close binary, and we include this high-order-multiple system 
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Fig. 2. — The properties of the target sample based on SDSS photometry (and thus not including the effects of close companions). Clockwise from top left: 1) 
the total mass of the binaries; 2) the mass ratios of the binaries; 3) the minimum resolvable separation (Rayleigh criterion for each telescope); and 4) the projected 
separation of the binary components. 
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Fig. 3. — A reduced proper motion plot for our targets. Primaries are blue 
and secondaries are red; the targets which are in fact high-order-multiples 
are plotted as triangles. Lines of tangential velocity for the disk, halo and an 
intermediate value are plotted for comparison. 

in our sample. Most of the data for SLW 1624-1-3249 was lost 
to an unfortunate disk problem. Data for that target presented 
here is on the basis of the few frames recovered, although 
much higher quality images were seen during observations. 

All observations were dark subtracted, flat-fielded, sky sub- 
tracted and aligned using a custom Python pipeline. The Palo- 
mar NGS-AO observations had very variable Strehl ratios due 
to the extreme faintness of our targets for natural guide star 



Fig. 4. — Achieved contrast ratios for the samples, for a companion detection 
5(T above both speckle and sky background noise. 

observations. To improve the image quality we took a few 
tens of 1 .4 second exposures for each target. Us ing the AO 
frame selection pipeline described in iLaw et akl C2009 ), we 
aligned each frame using the position of the brightest speckle 
in the point spread function (PSF), and selected the frames 
with the highest Strehl ratios to form a final co-added image 
with much improved resolution. 

3.1. Companion detection 
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TABLE 1 

The wide M-dwarf binary sample. 



Name 



SpTA 



SpTB Dist.[pc] Sep.[AU] Obs. Date Min. sep.[AU] Comp.sep.[AU] 



High-resolution H Keck LGS-AO 



SLW1558+0231 


15 58 51.5 


+02 31 14.7 


M3.0 


15 58 53.2 


+02 3136.9 


M3.5 


88 


3003 


2008/08/08 


2.4 


22±6 


SLW1602+08I2 


16 02 01.4 


+08 12 03.0 


M2.5 


16 02 01.0 


+08 11 55.6 


M4.0 


153 


1506 


2008/08/08 


6.1 




SLW1643+2100 


16 43 56.6 


+2100 37.8 


M2.0 


16 43 56.3 


+2100 31.3 


M4.0 


165 


1282 


2008/08/08 


6.6 




SLW1714+2807 


17 14 04.8 


+28 07 05.5 


M3.5 


17 14 04.2 


+28 07 04.1 


M5.0 


124 


909 


2008/08/08 


4.9 


10±2 


SLW2 1 37+0032 


21 37 05.4 


+00 32 31.3 


Ml.O 


21 37 04.6 


+00 32 46.1 


M3.5 


185 


3473 


2008/08/08 


7.4 


42±11 


SLW2I49+01I0 


2149 47.0 


+01 10 35.1 


M1.5 


2149 47.1 


+01 10 22.0 


M3.0 


173 


2279 


2008/08/08 


6.9 




SLW2I56+00I7 


2156 48.3 


+00 17 00.9 


M3.0 


2156 47.1 


+00 16 53.8 


M5.0 


143 


2807 


2008/08/08 


5.7 


27.2±6.8&11±3 


SLW23 19-0 105 


23 19 53.2 


-01 05 25.4 


M2.0 


23 19 52.3 


-0105 06.5 


M3.5 


173 


4056 


2008/08/08 


6.9 












Lower-resolution K^ Keck LGS-AO 










SLW0846+5407 


08 46 57.7 


+54 07 55.1 


M3.0 


08 46 58.7 


+54 07 46.8 


M4.0 


132 


1579 


2009/04/13 


7.2 




SLW0934+3425 


09 34 52.7 


+34 25 08.7 


M3.0 


09 34 53.3 


+34 25 14.4 


M3.5 


187 


1677 


2009/04/13 


10.3 




SLW1022+5835 


10 22 45.8 


+58 35 58.0 


M3.5 


10 22 45.1 


+58 35 50.8 


M4.0 


77 


700 


2009/04/13 


4.2 




SLW1053+0728 


10 53 09.0 


+07 28 35.2 


M2.5 


10 53 08.6 


+07 28 42.8 


M2.5 


101 


923 


2009/04/13 


5.6 




SLW1325+2253 


13 25 33.3 


+22 53 40.0 


M4.0 


13 25 32.4 


+22 53 57.9 


M4.0 


67 


1477 


2008/05/21 


3.7 




SLW1357+I952 


13 57 14.4 


+ 19 52 43.1 


M4.0 


13 57 13.1 


+ 19 52 22.5 


M4.5 


91 


2555 


2008/05/21 


5.0 




SLW1406+1546 


14 06 56.1 


+ 15 46 32.6 


M2.5 


14 06 54.8 


+ 15 46 05.8 


M3.5 


126 


4186 


2008/05/21 


6.9 


17±4 


SLW 1526+27 18 


15 26 48.3 


+27 18 29.8 


M4.5 


15 26 49.3 


+27 18 32.8 


M4.5 


58 


798 


2008/08/08 


3.2 




SLW1548+0443 


15 48 17.7 


+04 43 33.9 


M3.0 


15 48 18.3 


+04 43 39.4 


M3.5 


137 


1500 


2008/08/08 


7.6 


»400 AU (unconf.) 


SLW1632+I541 


16 32 13.3 


+ 15 4139.0 


M3.5 


16 32 12.5 


+ 15 42 03.8 


M3.5 


112 


3055 


2008/08/08 


6.2 




SLW 1639+4226 


16 39 50.6 


+42 26 12.1 


M3.0 


16 39 53.3 


+42 26 15.4 


M4.5 


148 


4426 


2008/08/08 


8.1 


120±30 


SLW2 1 27-0040 


2127 49.5 


-00 40 25.5 


M2.5 


2127 51.3 


-00 40 09.7 


M3.5 


203 


6342 


2008/08/08 


11.2 














Kp Palomar 200" NGS-AO 












SLW1049+6522 


10 49 08.8 


+65 22 25.2 


M4.0 


10 49 09.0 


+65 2157.4 


M3.5 


102 


2858 


2009/07/08 


11.3 




SLW1323+4429 


13 23 34.1 


+44 29 55.1 


M1.5 


13 23 31.9 


+44 29 49.2 


M2.5 


154 


3827 


2009/07/08 


17.0 




SLW 134 1-0222 


13 4102.5 


-02 22 00.3 


M3.5 


13 4102.8 


-02 21 48.2 


M4.5 


52 


662 


2009/07/08 


5.7 




SLW1345+2857 


13 45 26.0 


+28 57 23.7 


M2.5 


13 45 27.5 


+28 57 24.0 


M3.0 


106 


2159 


2009/07/08 


11.7 


12±3 


SLW1424+1856 


14 24 15.0 


+ 18 56 31.7 


M2.0 


14 24 14.1 


+ 18 56 25.6 


M2.5 


173 


2585 


2009/07/08 


19.0 


74±18 


SLW1427+0313 


14 27 22.6 


+03 13 49.4 


M3.5 


14 27 23.5 


+03 13 39.1 


M3.5 


85 


1452 


2009/07/08 


9.3 




SLW1439+5154 


14 39 32.5 


+5154 22.1 


M2.0 


14 39 33.8 


+5154 19.2 


M3.0 


122 


1590 


2009/07/08 


13.4 




SLW1500+2306 


15 00 47.5 


+23 06 26.9 


M2.0 


15 00 46.8 


+23 06 51.3 


M3.0 


140 


3684 


2009/07/08 


15.5 




SLW1512+2952 


15 12 07.1 


+29 52 05.8 


M3.5 


15 12 05.9 


+29 52 29.5 


M3.5 


60 


1710 


2009/07/08 


6.6 




SLW1521+5312 


15 2126.4 


+53 12 12.4 


M2.5 


15 2127.2 


+53 12 04.4 


M2.5 


88 


967 


2009/07/08 


9.7 




SLW 1522+5408 


15 22 21.2 


+54 08 27.9 


M1.5 


15 22 21.1 


+54 08 00.0 


M3.0 


120 


3343 


2009/07/08 


13.2 




SLW1530+3935 


15 30 21.9 


+39 35 22.1 


M2.0 


15 30 22.8 


+39 34 52.7 


M4.0 


149 


4628 


2009/07/08 


16.4 


26.2±6.6 


SLW1624+3249 


16 24 26.7 


+32 49 57.0 


Ml.O 


16 24 26.0 


+32 50 00.1 


M3.0 


133 


1373 


2009/07/08 


14.7 


64.2±16.1 


SLW1627+4336 


16 27 23.6 


+43 36 19.7 


Ml.O 


16 27 25.1 


+43 36 23.2 


M3.5 


160 


2582 


2009/07/08 


17.6 




SLW1352+4427 


13 52 18.1 


+44 27 09.7 


M3.0 


13 52 18.2 


+44 27 00.7 


M3.0 


103 


929 


2009/07/08 


11.3 




SLW 1546-0203 


15 46 41.5 


-02 03 06.4 


M4.0 


15 46 42.3 


-02 02 54.3 


M5.5 


55 


945 


2009/07/08 


6.0 





Note. — Positions are from the S PSS DR7 jAbazaiian et aH 2009); spectral types and distances are derived by fitting SDSS and 2MASS photometry using the 
spectral energy distributions given in Kraus & Hillenbran d ROOTa). Spectral types are accurate to approximately 1 subclass (extra precision is given to compare 
populations more easily). Photometric distances are accurate to approximately 0.3 magnitudes, or 15%. Min. sep. is the minimum detectable separation for an 
equal-m agnitud e close companion to either of the star in the target binary. SLW 1548+0443 A has a companion at a wide radius that requires CPM confirmation 
(Section l4!6!2l . 



Binary companions were searched for both automatically 
and by eye. In almost all cases the two components of the 
wide binaries were imaged simultaneously, giving excellent 
PSF references (within the limits of the isoplanatic patch). In 
other cases, the variable PSFs on these poor nights prevented 
exact matches with other targets in the survey, although com- 
parison with other stars often did allow rejection of persistent 
speckles. Three of the detected companions are at a contrast 
level where speckle noise is significant, and all those targets 
were imaged simultaneously with the other star in the wide 
binary. No candidate companions at those flux levels were 
seen in observations without direct PSF matches. 

The maximum contrast ratios for companion detection are 
shown in Figure |4] The photon and speckle background noise 
at a particular radius from the primary was determined by: 1) 
measuring the total flux in each of a set of circular apertures at 
that radius; 2) calculating the RMS variation in total flux; and 
3) requiring a 5cr companion detection above the calculated 
noise. 



3.2. Companion Photometry 

Aperture photometry was used to determine flux ratios be- 
tween the components of companions with separations >0.3". 
Closer companion flux ratios were derived by an iterative pro- 
cedure (following Kraus & Hillenbrand 2007b): 1) an esti- 
mate analytical point spread function is fit to each component 
of the binary; 2) the estimated PSF is used to subtract the flux 
from the secondary; 3) and a new PSF estimate is generated 
on the basis of the subtracted image. These steps are repeated 
until convergence, usually in 4-5 iterations. The final esti- 
mated PSF is fit to both components simultaneously to derive 
the flux ratio. Uncertainty estimates noted for the flux ratios 
include systematic uncertainties from this fitting process. 

4. RESULTS 

Ten of the 36 observed systems were found to contain close 
companions; one system is a quadruple consisting of two ap- 
proximately equal-brightness close binaries. The properties 
of the systems are summarized in Table |2] and images of the 
systems can be found in Figure |5] 
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The discovered companions 












Name 


Sep. ["] 


P.A. [deg] 


Dist. [pc] 


Sep. [AU] 


AK 


AH 


(ni/c)A/B 


SpTAa/Ba 


SpTAb/Bb 


Wide? 


SLW1345+2857Aa/b 


0.11±0.01 


51±8 


110±30 


12±3 


0.3 ±0.1 


0.31 ±0.05 


10.5 


M2.0 


M2.5 




SLW1406+1546Ba/b 


0.11±0.01 


79±8 


150±40 


17±4 


0.1 ±0.2 




12.1 


M3.5 


M4.0 


* 


SLW1424+1856Aa/b 


0.39±0.01 


276±2 


190±50 


74±18 


0.10 ±0.05 


-0.05 ± 0.05 


11.5 


M1.5 


M2.0 




SLWl 530+3945 Aa/b 


0.16±0.01 


94±5 


160±40 


26±7 


1.0 ±0.3 




11.5 


M1.5 


M3.5 


* 


SLW1558+0231Aa/b 


0.21±0.01 


229±4 


100±30 


22±6 


3.5 ± 0.3 


3.5 ± 0.3 


11.1 


M2.5 


M7.0 


* 


SLW1624+3249Ba/b 


0.38±0.05 


236±11 


170±40 


64±16 


0.7 ± 0.2 




12.0 


M2.5 


M4.0 




SLW1639+4226Aa/b 


0.59±0.01 


133.8±1.4 


200±50 


120±30 


0.63 ± 0.05 


0.69 ± 0.05 


12.4 


M3.0 


M3.5 


* 


SLW1714+2807Ba/b 


0.06±0.01 


259±14 


160±40 


9.7±2.4 


0.4 ± 0.2 


1.46 ±0.05 


13.5 


M4.5 


M5.5 


* 


SLW2137+0032Aa/b 


0.21±0.01 


122±14 


200±50 


42±11 


1.60 ±0.05 


1.69 ±0.05 


11.9 


Ml.O 


M5.0 




SLW2156+0017Aa/b 


0.13±0.01 


126±6 


210±50 


27±7 


0.4 ±0.1 


0.55 ± 0.05 


12.4 


M2.5 


M3.5 


^ 


SLW2156+0017Ba/b 


0.06±0.01 


149±13 


190±50 


11±3 


0.3 ±0.1 


0.27 ± 0.05 


13.7 


M4.5 


M5.0 


^ 



Note. — (ms:)A/B is the (unresolved) 2MASS K-band magnitude of the close binary system. SpTAa/Ba refers to the estimated spectral type of the primary 
component of each close-binary system; SpTAb/Bb is the estimated spectral type of the newly discovered companion. Contrast ratios are given in H and K, where 
available, and are relative to the nearest member of the target wide binary. The spectral types of the close-companion's primary, the close companion, and the 
other component of the wide binary are given with an accuracy of approximately one spectral subtype. Systems are noted as "wide" if they cross the Reid et al] 
1200 li) empirical mass/separation cutoff. The spectral types and distances have been refit using the resolved photometry, and hence diifer slightly from Table[T] 
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Fig. 5. — The newly discovered high-order multiple systems. Only the close-binary components are shown. All images are orientated with North up and East to 
the left and each image is 0.8 arcseconds square. Images are in H with the exception of SLW1530+3945 and SLW1406+1546, both displayed in K. 



4.1. Probability of physical association 

Our imaging dataset contains three possible companions at 
separations >3", found in a total survey area of 57600 arcsec-. 
All have SDSS and 2MASS colors inconsistent with being 
low mass stars at distances similar to our targets, and all are 
sufficiently bright to have been picked up during any of our 
observation runs. The sky density of possible companions 
in our fields is thus a maximum of a:9xl0"^ stars/arcsec^ 
(at 99% confidence). Our total companion search area at 
<3"radius is 2007 arcsec-; we thus expect 0.035 false as- 
sociations in our dataset, or equivalently a 3% probability 
of a single false positive association. This limit is correct 
for companions with AM<4 (those detectable in all our ob- 



servations). Our faintest candidate, SLWl 548 +0443 Ab, has 
a somewhat lower formal association probability because it 
would not have been detected in our Palomar imag ing, an d 
we leave its confirmation for CPM followup (Section r4.6.2b . 

4.2. Companion Statistics 

Eleven of the 71 observed stars are actually close bina- 
ries, leading to a raw binary fraction of 16^3% (with a bino- 
mial uncertainty estimation). However, our common-proper- 
motion wide binary selection imposes several biases on the 
binarity fraction of our target sample: 1) the selection of 
photometrically-clean, non-extended SDSS objects removes 
some binaries from the sample; 2) binaries are brighter than 
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single stars at the same colors in unresolved photometry, and 
so our magnitude-limited sample se lects binaries from a la rger 
space volume than the single stars dBurgasser et al.ll2003h : 3) 
those binaries found at larger distances are also likely to have 
decreased proper motion, and so may be removed from our 
target sample (Law et al. 2008); 4) the widely-separated com- 
ponents of high-order multiple systems containing close pairs 
can appear to be at different photometric distances in unre- 
solved photometry, removing them from the sample, and 5) 
the most distant (and therefore the physically widest) binaries 
will be de-selected by our 7"minimum-separation cut. 

4.2.1. Biases affecting proper-motion-selected close binaries in the 

SDSS 

We selected our sample from a subset of SDSS sources 
which are photometrically clean (no error flags), unsaturated, 
not flagged as galaxies, and also have a clean proper motion 
determination. Close, barely-resolved binaries in the SDSS 
may appear as extended sources and so be detected as galax- 
ies, and this could bias our selection away from systems with 
such components. Similarly, proper motions are only reliable 
if the nearest neighbouring object in the SDSS is >7"away 
(^unn et al. 2004, 2008); we thus remove objects with nearby 
resolved companions, and this may similarly bias us against 
systems with close companions. However, these biases only 
operate on systems which can be resolved as extended by 
SDSS. 

To determine the maximum detectable close compan- 
ion separation for our binaries we started with the < 
2.5"separat ion binaries in the Washington Double Star Cat- 
alog (Mas on et al.l 120011) . We then applied our photometric 
and proper motion selection cuts to those binaries which are 
present and unsaturated in the SDSS. We find that the SDSS 
star/galaxy separation does an excellent job of avoiding flag- 
ging extended binaries as galaxies - very few of the binaries 
are removed at this stage. However, many of the binaries 
are detected as doubles by SDSS, while additionally some 
sources have another unrelated object within 7"(although the 
removal of those targets does not bias our selection toward or 
away from binaries). For the purposes of this test we only 
remove the binaries that have another source within 3". This 
separation is chosen to include our test group of binaries at 
<2.5"separations (with some margin), but allows us to only 
flag sources affected by their own (de-blended) companions. 

The results of this selection are shown in Figure|6] All bina- 
ries with <LO"separation pass the cut (because they are seen 
as point sources by SDSS), and the higher-contrast binaries at 
wider separations also pass, because SDSS does not detect the 
faint companions. Our survey is therefore not significantly bi- 
ased against binaries up to 1.0"separation, or xilOOAU at our 
median target distance, and is sensitive to finding new low- 
mass companions at much larger rad ii from the target stars 
(the candidate brown dwarf in Section l4.6.2l for example). 

4.2.2. Monte-Carlo bias correction 

To generate correction factors for the other sample- 
selection biases, we generated an ensemble of simulated wide 
binary systems in the spectral type, proper motion and dis- 
tance range of our sample, and then applied our sample se- 
lection criteria to produce simulated target lists. This proce- 
dure is designed to correct both for the individual biases noted 
above and for any interaction between the various effects. 

We uniformly distributed 50,000 binary, triple and quadru- 
ple systems in a 300pc radius sphere. We then simulated 
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Fig. 6. — Determining the maximum binary separation which is not removed 
by SDSS quahty cuts. The open circles are binaries which fail either the pho- 
tometric or proper motion quality cuts due to their close companion; the filled 
circles are those detected as single stars with good-quality proper motions. 

SDSS photometry for each binary component using the spec- 
tral energy distributions (SEDs) given in Kraus & Hillenbrang 
(2007a). We included close companions in the generated un- 
resolved photometry, selected from the distributions described 
below, and included a 2% photometric uncertainty in the mag- 
nitudes. Proper motions for the binaries were taken from 
a Gaussian distribution with 25kms"' width, typical for our 
thin-disk targets (Section[2|i. 

To simulate our sample selection procedure we first ob- 
tained a best-fit spectral type and distance for the generated 
photometry of the binary's components; this step includes the 
effects of the added unresolved close companions on the in- 
ferred spectral types. Because our simulated photometry does 
not cover the range of SED variations inherent to real stellar 
populations, we add a 0.3-magnitude Gaussian-distributed un- 
certainty in the fitted distance modulus. Finally, we apply the 
distance, proper motion, separation, and magnitude selection 
criteria discussed in Section|2l 

We consider a variety of distributions for the spectral types 
of each component of the multiple systems, by averaging our 
results between the extremes of a flat distribution between Ml 
and M4, and simply setting the spectral type to our median 
spectral type of M2.5. We also average over a range of wide 
binary separations distributions, between all systems being at 
3000 AU, and a flat distribution between 1000 and 9000 AU. 
Each combination of distributions is also included so, for ex- 
ample, we consider systems with a flat distribution of wide bi- 
nary component spectral types, but with all close binary com- 
ponents having spectral types equal to their primaries. 

Finally, we count the number of selected simulated binary, 
triple and quadruple systems, and compare them to find the 
relative size of the "effective volume" from which each is 
taken. This volume directly sets the bias-correction factors for 
each type of system. The distribution of the close companion 
spectral types is the only population characteristic that has a 
>10% effect on the bias-correction factors, with an approxi- 
mately 2x change in effective volume between a flat early-M- 
star distribution and an all-equal distribution. 

We find that triple systems with one unresolved close binary 
component are suppressed in our sample by between 1 .4x and 
3.1x compared to wide binary systems. This is primarily be- 
cause a wide binary component that actually consists of two 
unresolved equal-brightness bodies appears to be closer to us 
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than the other component of the wide binary - and so is cut 
out by our aggressive Icr distance match requirement. In con- 
trast, the frequency of selected quadruple systems, not sub- 
ject to the above effect, is between 0.8x and 1.2x the binary 
fraction, depending on the distribution assumptions. In the 
high-order-multiple fractions given below, we include the full 
range of bias-correction factors in our formal uncertainty es- 
timates. 

4.3. Bias-corrected high order multiple fraction 

We discovered 10 high order multiple systems in our 36- 
target sample, for a raw fraction of 28^g%. 

The bias correction factor for triple and quad systems dif- 
fers by almost 2x, and the detection of a triple system does 
not exclude the possibility that it also contains a further un- 
resolved fourth component. However, only one of the 10 
high-order-multiple systems is a detected quadruple system. 
Since quadruple systems have an approximately 2x greater 
chance of surviving our sample selection than triple systems, 
we can interpret this as evidence that quadruple systems are 
quite rare, at a <5% true incidence for our separation range. 

There is an alternative possibility: that a larger fraction of 
the wide binary systems are in fact quadruple systems, with a 
small probability that one or both close companions would be 
at wide enough radius for detection in our AO survey. Some 
indication of this may be found in the high-resolution sub- 
sample, where the only quadruple system in our sample is 
found with both components at small separations, but a much 
larger sample is required to fully evaluate this hypothesis. For 
the purposes of the rest of this paper we take the simplest 
approach: since 90% of the detected systems appear to be 
triples, we assume that the true population of multiple sys- 
tems is dominated by triple systems. This assumption changes 
only the absolute value of the high-order-multiple fraction, 
and does not affect detailed comparisons of multiplicity rates 
within our sample population. 

Adopting the triple bias correction factor from Section l4.2.2l 
for all 10 multiple systems, and using the binomial distribu- 
tion probability formalism detailed in lBurpasser et al.l (l2003h . 
the corrected high-order-multiple fraction is 45^jg%. 

4.4. Narrow-angle incompleteness 

We designed the high-resolution Keck H subsample to have 
sufficient resolution (< 4 AU at lOOpc) to probe the 4 AU 
peak of the field-binary separation distribution of mid-M- 
dwarfs. Three triple systems and one quadruple system were 
detected in this 8-system high-resolution subsample, leading 
to a high-order-multiple fraction of 70^27*^' statistically com- 
patible with the binary fraction derived in the lower-resolution 
portion of the survey. Indeed, only one (20%) of the com- 
panions detected in the Keck high-resolution survey would be 
undetectable in the lower resolution portion of the surveys. 
Although these are small number statistics, it appears that in- 
sensitivity to companions in the 4-10 AU range does not lead 
to a large change in the number of detected companions. 

4.5. Contrast-ratio incompleteness 

Eight of the eleven close companions are well above our de- 
tection sensitivity (Figure |8]l, and they are evenly distributed 
in the detectable separation range. The remaining three com- 
panions are detected just above the sensitivity limit. Although 
these are clear detections, we cannot exclude the presence of 
a large population of higher-contrast binaries. 




Fig. 7. — An //-band image of the faint possible companion around 
SLW1548+0443A. 
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Fig. 8. — The contrast ratios of the detected companions compared to our 
detection sensitivities. 

4.6. Detected binaries with a high mass ratio 
4.6.1. A very low mass companion to SLWl 558+0231 A 

All but one of the confirmed close binaries is approximately 
equal-mass. The exception, SLW1558-i-0231Ab, has contrast 
ratios of A// = 3.4 + 0.25 and A/T = 3 .5 ±0.25 and a pr o jected 
separation of 22.3+5.6 AU. Using the Bar affe et all ( Il998b 
models, and our calculated distance modulus for the system, 
at 5.0 GYr age the secondary is best modelled with a mass of 
0.090±0.007Mo (Icr uncertainties). At a 500 MYr assumed 
age the secondary mass is lowered to 0. 080+0. OO9M0. 

4.6.2. An unconfirmed wide brown dwarf around SLWl 548+0443 A 

A faint object is found just under 3" (~400AU) away from 
SLW1548+0443A (FigureHi. The object is ~7 mags fainter 
than the target M3 star and is not visible in SDSS because 
of the high contrast. If the object is associated with the 
wide M-dwarf binary it has absolute magnitudes My«14.9 
and M^sil3.7, which is consistent with a mid-mass brown 
dwarf at the target star distance. Although the probability 
of association is fairly high, CP M confirmation is required 
for this companion (Section [4. II ). The system's proper mo- 
tion is s;0.06"/yr ', and so it should be possible to con- 
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Fig. 9. — A histogram of the separations of the newly-discovered close com- 
panions; Icr Poisson eiTor bars are shown. 

firm the object's association on a short timescale. Since this 
candidate remains unconfirmed, we do not include it in any 
subsequent analysis. We will address its association with 
SLW1548-f0443A in more detail in a future paper. 

5. DISCUSSION 

5.1. The Population Statistics of Close Companions in Wide 
Binary Systems 



The high-order-multiplicity fraction of our sample, 45^[g%, 
is statistically compatible with the 25% fraction found in 
the narrower M -dwarf multiple systems in the 8pc sam- 
ple JReid & Haw lev 2005). Similarly, the presence of only 
one detected quadruple system is compatible with the single 
iReid & Haw lev (2005) quadruple detection. 

Our detected companions are almost all at low contrast ra- 
tios and relatively small separations on the sky, although two 
possible very low mass companions were found at larger radii. 
The distribution of mass ratios is compatible with that found 
among field close bi naries in this mass range (Reid et al] 
l200ltlLawetani2008l) . 

The close-binary component separation distribution is 
strongly peaked at separations <30AU (Figure|9]l; this is three 
times closer than our ~100AU maximum separation cutoff 
(see Section |4. 2. lb . At <30 AU separations too few systems 
are detected to allow detailed analysis, but it is worth not- 
ing that we do not see see a strong peak towards even closer 
systems. The survey is sensitive to companions at separa- 
tions <20AU around all our targets, and <10AU around two 
thirds of our targets. A very strong peak towards closer sepa- 
rations, such as that observed in lower-mass M-dwarfs in the 
field (Law et al.,2006.i2008i:iSiegler et al..2005i) . would be de- 
tectable. 

5.2. High-order Multiplicity 

Next, we assess if the properties of the target wide binaries 
have an effect on the probability of high-order-multiplicity. In 
the following discussion, we refer to the newly detected close 
companions as "close binaries", and the systems of which 
they are part as "wide binaries". Binning the targets by the 
separation of the wide binaries, we find a marginally signifi- 
cant (2cr) increased high-order-multiple fraction for the wider 
targets in the survey (FigurefTOb. For systems with separations 
up to 2000AU the high-order-multiple fraction is 21+^^%, 
compared to ll^j^^'^c for systems with separations >4000AU. 



^ 






-a 
o 



60 



80 














60 












40 










20 















2000 4000 

Wide binary separation / AU 



6000 



Fig. 10. — High-order multiplicity as a function of the separation of the wide 
binary systems. The high-order- multip le fractions, bias corrections, and error 
bars are calculated as in Section 14.21 with appropriate Monte-Carlo separa- 
tion distributions for each of the bins. To compare the different separation 
ranges without including the range of bias coirections' extra uncertainty in 
the absolute high-order-multiple f racti ons, we average over all the bias cor- 
rection models detailed in Section|4J2] The final bias corrections at the three 
separation ranges differ only at the 5% level. 
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Fig. 11. — The median close companion separation as a function of wide 
binary separation. We caution that the leftmost bin contains only two data- 
points (10 AU and 64 AU) and should not be interpreted as evidence for a 
varying median separation. 

Indeed, three of the four widest binary targets in the survey 
are high-order-multiples, while we find only one triple sys- 
tem from eight <1000AU systems. 

There is no evidence that the median close binary sepa- 
ration changes as a function of the wide binary separation 
(Figure [TT]). This suggests that the variations in high-order- 
multiplicity are not an effect of detection efficiency or tar- 
get selection changing with system distance (as would be ex- 
pected if uncorrected biases remain in our binarity statistics), 
but rather are a real change in the multiplicity fraction. 

Of the 23 unusually wide systems, six have close-binary 
components, com pared to fou r close-binary components in 13 
systems below the iReid et all ([2001) cutoff (Figure[T2]i. With 
bias-corrected high-order-multiplicity fractions of 42.7'^^}j'\% 

and 49.0;^2o 7*^ respectively, the close-binary fractions are sta- 
tistically indistinguishable. The close-companion mass ratios 
of the two samples are similarly indistinguishable, while the 
22 + 7 AU median separation of the unusually wide binaries' 
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Fig. 12. — The separations and masses of our wide binary sample {circles); 
our detected high-order multiple systems {filled circles) and the young wide 
binary sample from Kraus & Hillenbrand (2009) {triangles). The effect of 
binarity on the estimated total system mass is shown for each of our detected 
higher-order multiple systems, with the lower points based on unresolved 
photometry and the upper points on AO resolved ph otometry. The empirical 
field binary separation cutoff from lReid et alj 120011) is shown as a solid line. 
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Fig. 13. — A comparison of the estimated binding energies for our systems 
before and after the detection of extra components. Initial binding energies 
from unresolved photometiy are shown for all our targets; an'ows point to 
the binding energies after finding additional components. The dashed fine is 
designed to bisect our sample, and h as the same binding -energy/mass slope 
as the Jeans-length-based relations in lFahertv et aU2010[ 

close-companion population is not significantly different from 
the 53 + 24 AU normal-binary median separation. 

Figure [13] shows the effect of the resolved companions 
on the systems' binding energies. Most of our detected 
high-order-multiple systems are found towards the bottom of 
the figure, at low binding ener gies and masses. A lthough 
all these systems are above the ' Fahertv et al.l ( 1201 Oh Jeans- 
length-based binding energy cutoff, it is instructive to con- 
sider if there is evidence for the wide binary system proper- 
ties changing in a similar form to that relation. In Figure [T3] 
we also plot a curve with the same power-law slope as the 
iFaherty et al.l (12010.) relation, with an additional multiplica- 
tive factor so that the binding energy cutoff bisects our sam- 
ple. The full form we adopt is E = 20M^'^, where E is the 
binding energy in units of lO'*' ergs, and M is the system total 
mass in solar masses. There are two high-order-multiple sys- 
tems above this line (a bias-corrected 20^g^%), and eight be- 
low the line (62^\\%). Although we must be cautious to avoid 
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Fig. 14. — The frequency of high-order multiple systems as a function of 
total system mass. The system mass of the close binary systems is estimated 
here from uni'esolved photometry, for direct comparison with the systems 
without detected close companions. There is no significant evolution of the 
high-order-multiple frequency across this mass range. 

simply picking a relation which fits the data, this marginally 
signficant difference, and the lack of a similarly significant 
difference for a simple binding-energy or mass cutoff (Fig- 
ures [13] & [T4l i. may suggest that Jeans-length considerations 
are indeed important in the formation of high-order-multiple 
wide M-dwarf systems. Much larger sample sizes and more 
detailed statistics are required to fully investigate this. 

5.3. Formation scenarios 

Our results demonstrate that the individual components 
of unusually wide M-dwarf binary pairs are broadly simi- 
lar to other low-mass stars in the field, as both populations 
have similar companion frequencies, separation distributions , 
and mass ratio distributions (e.g. iFischer & Marcyl Il992t 
iSieder et aUlIOOllLaw et ani2008h . 

Less than half of the systems in our observed sample have 
additional companions at a separation of >5-10 AU, a com- 
parable rate to typical field M dwarfs. We are not sensitive to 
companions at separations of <5-10 AU, so it is possible that 
the systems have unrecognized additional components inside 
the detection limits of our survey. However, the absence of a 
significant additional population at separations of ~5-10 AU 
in our high-resolution subsample seems to argue against this 
hypothesis. 

The fact that many of our observed systems appear not to 
contain extra companions (at our detection limits) strongly 
suggests that an unusually wide separation is not necessar- 
ily an indication of (unrecognized) excess mass; some wide 
binary systems appear to genuinely violate the a^ax - Mtot re- 
lation seen for most low-mass stars. As this relation was set 
empirically on the basis of the then-discovered population of 
M-dwarf binaries, this is not in itself surprising. 

However, we do find a 2cr increase in the frequency of high- 
order multiplicity for the widest pairs. This is an intriguing 
result, as both our survey's smallest-separation and largest- 
separation binaries would be considered unusually wide com- 
pared to the vast majority of field binaries (Reid et al. 2001), 
and no current models predict that there should be such an 
evolution in wide binary properties. One possible explana- 
tion is that the very widest binaries do in fact need higher 
total masses to form or to survive, leading to a higher proba- 
bility of finding extra components in the systems, although it 
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remains to be explained why there would be a rapid change 
in binary properties at such large separation ranges. Larger 
sample sizes are required to confirm this observation. 

Surveys of star-forming regions show that the frequency 
of wide (>500 AU ) binary pairs declines smoothly with 
system mass ( Kraus & Hillen brand 2009), and we find no 
evidence that the high order multiplicity of our systems 
changes with total system mass. It therefore seems plau- 
sible that perhaps all protostellar cores undergoing freefall 
collapse could fragment into multiple systems (as reviewed 
by Bodenheimer & Burkert 2001) on any length scale, but 
the probability of this collapse occurring at a given length 
scale (i.e. the probability of achieving a Jeans-critical over- 
density) declines with decreasing characteristic density, and 
hence with decreasing total mass of the core and the resulting 
binary system. The two resulting fragments would be oth- 
erwise unremarkable, so they would then evolve and further 
fragment in a manner like low-mass protostars that are not in 
wide binary pairs. 

This idea would explain both the rarity of wide systems 
and the similarity between their components and other low- 
mass stars. However, other proposed models could also match 
the observations, such as the mutual ejection hypothesis of 
iKouwenhoven et aH (12010 ). which postulates that wide binary 
systems could form via simultaneous ejection from the natal 
cluster on a similar trajectory. The validity of such a mech- 
anism should be tested as part of large-scale star formation 
simulations (e.g. Bate 2009; Offner et al. 2009). The increase 
in high-order-multiplicity for the very widest binaries, if con- 
firmed, could offer a sensitive new test for formation models. 

6. SUMMARY 

We find 11 new close companions in our survey of 36 
M-dwarf binaries with separations >500AU, giving a bias- 
corrected high-order- multiple fraction of 45^ [g% for these 
wide systems. The close-binary fractions and close-binary 
mass ratios of our wide binary targets are statisticall y indis- 
tinguishable whether they are above or below the Rei d et all 
(12001) empirical "unusually wide" mass cutoff. We do, 
however, find marginally significant evidence for an evo- 
lution of high-order-multiplicity following the form of the 



iFahertv et alj ( 120101) Jeans-length binding energy relations. 

The unusually wide binaries' close-companion separations 
are not significantly different from the normal-binary sam- 
ple. The close-binary component separation distribution is 
strongly peaked at separations <30 AU, but we do not see see 
a strong peak towards even closer systems. Almost all the 
detected companions have similar masses to their primaries, 
although two very low mass companions, including a candi- 
date brown dwarf, were found at relatively large separations. 

We find 2cr evidence for an increased high-order-multiple 
fraction for the widest targets in our survey, with a high-order- 
multiple fraction oi21*]^% for systems with separations up to 

2000 AU, compared to 77^22"^ fo'" systems with separations 
>4000AU. Three of the four widest binary targets in the sur- 
vey are high-order-multiples, while we find only one triple 
system from eight systems with a wide separation <1000AU. 
It is possible that the widest binaries do in fact need higher 
total masses to survive, leading to a higher probability of find- 
ing extra components in the systems. Larger sample sizes are 
required to confirm this. 

Our results indicate several potentially exciting paths for- 
ward in studying unusually wide binary systems. The high 
frequency of high-order multiples of the widest systems in 
our sample should be confirmed with larger surveys, and even 
wider systems should be surveyed to see if the frequency con- 
tinues to increase. Our sample will also be pushed to much 
lower total system masses under better LGS observational 
conditions, and the existing sample should be surveyed with 
high-dispersion spectroscopy to search for extremely close bi- 
nary companions. 
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